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Abstract 

<N 

5_i We study a class of noncommutative gauge theory models on 2-dimensional Moyal 

Qh space from the viewpoint of matrix models and explore some related properties. Ex- 

panding the action around symmetric vacua generates non local matrix models with 

polynomial interaction terms. For a particular vacuum, we can invert the kinetic oper- 
^vq ator which is related to a Jacobi operator. The resulting propagator can be expressed 

in terms of Chebyschev polynomials of second kind. We show that non vanishing cor- 
' _| ' relations exist at large separations. General considerations on the kinetic operators 

-4— > stemming from the other class of symmetric vacua, indicate that only one class of 

symmetric vacua should lead to fast decaying propagators. The quantum stability of 
1) the vacuum is briefly discussed. 
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1 Introduction. 



In Noncommutative Geometry (NCG) [I], one of the guiding ideas is to generalize the 
duality, existing in Riemannian geometry, between spaces and associative algebras in a 
way that the structural properties of the space, for instance topological, metric, differential 
properties, can be given an algebraic description. It turns out that many of the building 
blocks of modern physics fit well with the basic concepts of NCG which may ultimately lead 
to a better understanding of spacetime at short distance. For instance, NCG provides a 
way to resolve the physical objections, emerging from the concurrency of General Relativity 
and Quantum Mechanics [2], to the existence of continuous space-time at Planck scale. 
Once the noncommutative nature of space(-time) is assumed, it is natural to consider 
field theories on noncommutative spaces, called Noncommutative Field Theories (NCFT). 
NCFT emerged around 1986 in String field theory [8], followed a few years later by the first 
NCFT's on the fuzzy sphere, a finite dimensional NCG, [3]. NCFT on the Moyal space 
were singled out as effective regimes of String theory 0[B] around 1998 and as underlying 
structures in quantum Hall physics [Tj, |8]. For reviews on Moyal NCFT, see e.g [9]. 

The renormalization of NCFT is difficult, unless the underlying NCG is finite. In Moyal 
geometry, this is due to the ultraviolet/infrared (UV/IR) mixing, appearing already in 
the real- valued ip 4 model on the 4-d Moyal space [TO]. It comes from UV finite nonplanar 
diagrams which exhibit IR singularities generating new divergences when inserted into 
higher order diagrams that cannot be cured. A first solution to the mixing in the scalar 
field theory amounts to add to the initial action a harmonic oscillator term. This yields 
the Grosse-Wulkenhaar model, perturbatively renormalizable to all orders |11| . |12j . [13] . 
Various related properties have been studied, among which classical and/or geometrical 
ones, 2-d fermionic extensions [2]-p2]. The Grosse-Wulkenhaar model has vanishing of 
the /3-function to all orders |18j when it is self-dual under the Langmann-Szabo duality 
[19 1, and is very likely to be non-perturbatively solvable [201. Scalar field theories on the 
noncommutative space R?, a deformation of M 3 , which are free of UV/IR mixing have 
been built recently in (21] . Whether this simply comes from the low "dimension" of the 
space or reflects a specific property of the underlying NCG remains to be seerQ . 

The UV/IR mixing also occurs in gauge models on 4-d Moyal space [23] . For early 
studies, see e.g [24J and references therein. The mixing appears in the naive noncom- 
mutative version of the Yang-Mills action given by So = \ f d 4 x(F^ u ★ F lll/ ){x), showing 
up at one-loop order as a hard IR transverse singularity in the vacuum polarization ten- 
sor. Attempts to extend the Grosse-Wulkenhaar harmonic solution to a gauge theoretic 
framework have singled out [26] a gauge invariant action expressed a^] 

d 4 xl K -F llu *F llu + —{A ll ,A u }l + KA fl *A fl ). (1.1) 

where f2 and k are real parameters, while A^ = + ^x^ is a gauge covariant one- form 
given by the difference of the gauge connection and the natural gauge invariant connection 
induced by the minimal (e.g. based on the derivations d^) Moyal differential calculus [25] 
(see section]!]). The corresponding mathematical framework has been developed in [27] 128]. 



Eq. (1.1) can also be viewed as a spectral action related to a so called finite volume 
spectral triple [29], whose Dirac operator is a square root of the kinetic operator of the 
4-dimensional Grosse-Wulkenhaar model. The related noncommutative (metric) geometry 
is rigidly linked to the Moyal (metric) geometries, as shown in |30 [ I31 ] |32]. Unfortunately, 



X A discussion on the origin of the mixing for translation-invariant products can be found in [22] 
2 Notations and conventions are collected in the section [2] 
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the action ( 1.1 ) is hard to deal with when it is viewed as a functional of the gauge potential 

, Sfi [Afj] . This is mainly due to its complicated vacuum structure explored in [33] which 
excludes the use of any standard perturbative treatment. Other attempts to control the 
UV/IR mixing have been considered in [33] -[38]. However, showing that any of these 
models is renormalizable is still an open problem. 

When expressed as a functional of the covariant one- form A^, the action (1.1) bears 
some similarity with a matrix model, where the field An can be represented as an infinite 
matrix in the Moyal matrix base (see Appendix [B]) . To our knowledge this interpretation 
has not been explored so far for the action (1.1), although the matrix model formulation 
of NC gauge theory has been known since many years |39j . Keeping in mind that A^ is a 
natural variable of the Moyal geometry [27] , a sensible question is to explore the properties 
of the action as a functional of the field An , in order to determine to what extent Sq [An] 
may give rise to a meaningful quantum theory. 

This requires at least to choose a vacuum, expand the action around it, and perform 
the difficult but mandatory computation of the propagator. In the 4-d Moyal case, an 
additional complication comes from the need to control the UV/IR mixing behavior of the 
ghost sector coupled to the gauge sector. This can be easily overcome in the 2-d Moyal 
case, that we will consider in this paper, thanks to a suitable gauge choice, akin to the 
temporal gauge, for which the ghosts decouple from the gauge sector. 

The purpose of this paper is to perform a first exploration of the gauge theory on 
2-dimensional Moyal space described by the action of Eq. ( 1.1 ), at Q ^ 0, viewed 

as a matrix model as explained above. The expansion of the action around the non-trivial 
symmetric vacua determined in [33] , to which we restrict in this paper, and a BRST further 
gauge-fixing give rise generically to a non local matrix model with a complicated kinetic 
operator together with cubic and quartic polynomial interaction terms. For a particular 
symmetric vacuum corresponding to Q = |, we show that the kinetic operator is a Jacobi 
operator. The computation of the propagator can then be carried out. The resulting 
propagator can be expressed in terms of Chebyschev polynomials of second kind. We 
show that non vanishing correlations exist at large separations. The quantum stability of 
the vacuum is briefly discussed in the particular situation. A divergent 1-point function 
appears signaling likely a loss of symmetry at the quantum level. 

Extending our analysis to general considerations on the kinetic operators stemming 
from the other vacua found in [33], we single out a particular class of symmetric vacua 
that should lead to fast decaying propagators at large separation, reminiscent of a quasi- 
local behaviour for the corresponding matrix models. 

In section[2]we collect relevant properties of Moyal geometry that we will need together 
with specific features of the action ( 1.1 ) at fi = 0, such as the absence of UV/IR mixing and 
additional cancellation between IR singularities specific to 2 dimensions. The case f2 ^ 
is considered in section [3] which involves the general expansion and gauge fixing together 
with the computation of the propagator. We finally discuss the results and conclude. 



2 Yang-Mills type theory on R^: The Q = case. 

2.1 Basics on Moyal algebra and gauge invariant actions. 

In this subsection we collect useful properties of the Moyal algebra. For more details, see 
e.g [ID]. Let 5(M 2 ) = S and 5'(M 2 ) = S' be respectively the space of complex- valued 
Schwartz functions and the dual space of tempered distributions on IR 2 . The associative 
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Moyal ^-product is defined for all /, g in S by the map: * : S x S — > S 

(/ * 5)(x) = PF / d2yd2z f^ + y)9{x + z)e-^ e ^ z \ 0^ = 0^ J) (2.1) 

where 6 G R, 9 > The integral is a faithful trace: J d 2 x (f * g)(x) = f d 2 x (g* f)(x) = 
f d 2 x f(x)g(x). The Leibniz rule holds: d^{f * g) = d^f * g + f * d^g, V/, g G S. The 
complex conjugation in 5, a 4 a', Va 6 5, defines an involution in S. It extends to 5' 
by using duality of vector spaces. The ^-product (2.1) can be extended to S' x S — > S' 
using again duality ((T * a,b) = (T, a*b), VT G 5', Va, b G S) and continuity of (2.1 ). In 
a similar way, (2.1 ) extends to S x S' — > S', via (a *T,b) = (T, b ★ a), VT G 5', Va, b & S. 
Here, the Moyal algebra denoted by Eg is the multiplier algebra of (5,*), := A^Pl-M/? 
where M i = {Te57fl*Te5, Va G 5} and Mr = {T G S' / T * a G S, Va G 5} are 
respectively the left and right multiplier algebras of (5,*) |40j. We set [a, b]+ := a-kb — b-ka. 
For any a, b G M^, the following relations hold true: 

d^(a*b) = d^a-kb + a-k d^b, {a * b)^ = * a) , (2.2) 

x^-ka = (x^ ■ a) + -Q^vdud, a -k x^ = (x M • a) - -@f, u d u a. (2.3) 



Note that Eq. (2.3) implies the celebrated relation [x^, x„]* = iQ^ among the coordinate 
functions of Rg. In section [5J we will use the matrix base of M. 2 , [30]. The relevant material 
is given in appendix [B| 

As shown in |27| , the algebraic properties of classical gauge invariant actions on the 
Moyal plane are described by a simple version of the derivation-based differential calculus. 
Here, the unital algebra itself , M.g, plays the role of the right Mg-module while the Lie 
algebra of derivations is the Abelian algebra generated by the derivatives d^, /i = 1,2. 
The map Vx : ^ -> Kg, X = (fy), /i = l,2 given by 

= V M (a) =d fM a-iA fl *a, ; = iV M (l K 2), Va G (2.4) 

with A* = defines a hermitian connection, V, for the hermitian structure h : Kg xIR^ — > 
Kg, h{mi,m2) = m\ *m,2,Vmi,m2 G Ra. The curvature of V is 

R(X, Y) : R 2 g ^ Rq, R(X, Y)m = ([V x , Vy] - V[x,y])m, Vm G Eg (2.5) 

for any derivations X and Y and in the present case yields 

R„» ■= R(d», d v ) = -i(8^A u - d v A^ - i[A^ A u ]+). (2.6) 

To make contact with usual notation, we set from now on R^ v = —iF^ v . The unitary 
gauge group U(M.g) acts as V 9 X = g^ o Vx 9- This yields 

Al = gUA^g + igUd^g, F* =g**F lu ,*g, VgeU{R 2 e ). (2.7) 



We define x M := 2Q^x u . Recall [27] that the one-form components A^ = — \xfx := A^ 
define a gauge invariant connection one-form so that the map (2.4) becomes: 

ViT(a) = 3 M a + * a = |a* x M , Va G Kg, (2. 
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We will use the notation yO z s y M M „ z" and Einstein summation convention. 
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where gauge invariance (V 1 ™) 9 = gn o Vy" o g = V^ 11 follows from the second equality in 
Eq. (2.8). The corresponding curvature is F % ™ = —iQ^. Thanks to the existence of the 



gauge invariant connection, a natural covariant one-form can be defined as 



A„ := i(V„ - V;T) = A„ + -£ MJ ^ = ^*^*^ V 5 eW(Kg). (2.9) 



.A^ are formally the "covariant coordinates" introduced in [JT] on Moyal space-time in the 
presence of a local symmetry, although here they play a different role, namely that of the 



dynamical gauge variables. From Eq. (2.9) and the first of Eqs. (2.3), one obtains 



F^ u = -i[Ap,A v ]i, + e;J. (2.io) 

The extension to nSNis straightforward. For n = 2, the above framework underlies 



the action (1.1) which can be expressed solely in terms of A^ |26j, namely: 

S n = [d 4 x [-^[A^A^l + ^iA^A^l + KA^A^. (2.11) 
The gauge invariance is obvious in view of the tracial property of the integral and the 2nd 



relation in (2.9). Note that this action can be obtained from a spectral triple of a specific 
type introduced in [29] and shown, in the sense of noncommutative metric spaces, to be 
homothetic [30] to the standard Moyal spectral triple. 

2.2 The Q, = case. 



When the real parameters Q and k are set to zero, the action (2.11), expressed as a 
functional of An, exhibits UV/IR mixing which shows up at the one- loop level as a hard 
IR singularity in the vacuum polarization tensor. The situation is more favorable for the 



2-d version of (2.11). The corresponding model on at fi = k = is known to be 
UV/IR mixing free. Only the sector of planar diagrams actually matters and this latter 
is similar to the perturbative expansion of a commutative 2-d Yang-Mills theory. In fact, 
by using a "temporal-like" gauge, e.g A2 = 0, the gauge fixed action is purely quadratic 
and separates into a gauge potential A^ part and a ghost part, as for commutative 2-d 
Yang-Mills theories. Alternatively, one may use the popular "covariant" Landau gauge 
for which the gauge fixed theory "does not look free". By standard calculation, it can 
be easily realized that the hard as well as logarithmic IR singularities in the vacuum 
polarization tensor responsible for the UV/IR mixing are proportional to (2 — d) [27], [36] . 
This cancellation propagates to other higher order correlation functions as a consequence 
of Slavnov- Taylor identities [37J. We close this section by recalling that massless 2-d 
theories are known to exhibit additional IR singularities. In the case of 2-d commutative 
Yang-Mills theories, these depend of the gauge function. It is instructive to examine more 
closely the fate of these additional singularities within Sq=q L4J by using the Landau 
gauge. The gauge fixed action is Stat = <Sn=o [-^/J + Sgf with 

S GF = s J d 2 x(Cd»A^) = J d 2 x{bd»A^ - Cd^{d^C - i[A^ C]*)) (2.12) 

where the Slavnov operation s is defined by sA^ = d^C — iLA^, C]*, sC = iC*C, sC = 
b, sb = 0. Here, C, C and b are respectively the ghost, the antighost and the Stuckelberg 
field with ghost number equal to +1, —1 and 0. s acts as a graded derivation with grading 
defined by the sum of the degree of forms and ghost number (modulo 2) and s 2 = 0. 
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Standard calculation using (A. 4) and (A. 5) leads to the planar (uj p (p)) and non planar 



(uj^ p (p)) 1-loop contributions to the ghost 2-point function: 



(p) 



(2jr) D / 2 (p 2 



r 



r(f)r(f 



D-l) )L[Z 2 h 



(2.13) 



< P (p) = -2^/ q dxx (27T)D/ 2 



^(M 2 )f- 2 (^(v^) 2 ^^ 2 _ f (v / ^)), (2-14) 



For D = 2, the UV finite planar contribution (2.13) has an 



where M 2 := p 2 x(l — 
extra(2-d) IR singularity coming from the factor r(y — 1). Now, the UV finite non planar 

= (M 2 )- 1 , while the 
\ + (az + bz 3 + ...) + 



contribution ( 2.14| ) has an IR singularity coming from (M 2 



z) ~ 



factor (^p 2 M 2 ) 2 -^K 2 ^d{^p 2 M 2 ) is IR finite. Using K\ 

log(z) + ... where a, b G R, one checks that only the first term combines with (M 2 



to give an IR singularity. Therefore, the small \p\ behaviour of (2.14) is 

1 



1 

2^ 



o 



x(l — x] 



+ IR regular, p ~ 0, 



(2.15) 



This IR singular term is exactly cancelled by the planar contribution (2.13) at D = 2: 

■i 



w P [p) = — / dxx— r 

9 W 2vr J x(l - x) 



Hence the small p limit of cu g (p) is finite: 



limoj g (p) = ]im(u) g(p) +ui g yp (p)) = finite. 



(2.16) 



(2.17) 



p-»0 p— >o 

Similarly, the planar and non planar 1-loop contributions to the polarization tensor uo^ u (p) 



arc 



1 



dx 



1 



7T Jq p 2 x(l — X) 



(2.18) 



d D k ^(PA*) ((2 _ D)[ ^ 2fc + ^ (D + 2) (2 _ 1Q) 



W ^ W ~ "J (2n) D k 2 (k + p) 2 
Setting 0£> = 2 and using (IA.5I), we obtain 



1 



<W — PfiPv)- (2.20) 



7T 7 p 2 x(l — x) 

As above, by using the asymptotics given in the appendix [Al it can be checked that the 



extra IR singularity in (2.20) exactly cancels the planar contribution uj p v (2.18). Hence 



<*V(p) = 7r{p ){p 5^ - p^Pu), lim 7r(p ) = finite 

p->0 



(2.21) 



From (2.17) and (2.21), one concludes that IR singularities of 2-dimensional origin occur- 



ring in the 1-loop planar parts of ghost and gauge 2-point functions are exactly compen- 
sated by their non planar counterparts. This cancellation extends to the 3- and 4-point 



4 From ||A.6||, the 1st two terms between brackets in (2.191 behave respectively as log \p\ and 



f.1 I 



signaling UV/IR mixing unless D = 2. 
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functions as can be seen by using the Slavnov- Taylor identities [37]. On general grounds, 
one can expect that these 2-d IR singularities depend on the gauge choice. This can be 
easily exemplified by considering the case of 1-loop 2-point functions for a commutative 
2-d pure Yang-Mills theory, for instance either in the Landau gauge or in the temporal 
gauge, whose behavior is the same as the planar part of the above 2-point functions: IR 
singularities appear in the first gauge while they are simply absent in the other gauge. 
Notice that the present situation is slightly different since there are no such IR singulari- 
ties for both gauge choice. It suggests that the absence of these 2-d IR singularities does 
not depend on the gauge choice. The coupling to a fermion has been discussed in [37] and 
gives rise, as expected, to a dynamical mass generation mechanism as in the Schwinger 
model. 



3 Noncommutative gauge theory as a matrix model 



Many attempts to deal with the induced gauge theory on Moyal space have tried to 
interpret it as a Yang-Mills type theory. Indeed, the action was assumed to depend on the 
gauge potential or in other words the action (3.2) is considered as a functional of the 
gauge potential : S = S[A^]. Promoting this interpretation beyond the classical level is 
still unsolved. In this section, we will change the above viewpoint and use the covariant 
field Afj, defined in Eq. (2.9) as the fundamental variable entering the action. Doing this, 
we will therefore formulate the induced gauge theory (2.11) as a matrix model S = S[A^\ 
and examine if such a matrix model invariant under An — > Afi = g^ * A^ -kg, Vg G U(M>g) 
can have a consistent interpretation beyond the classical order. To our knowledge this 
interpretation has not been explored so far, despite the formal similarity between the first 



term of the action (2.11) and the (bosonic part of) the action for some type IIB matrix 
models, such as the IKKT matrix model considered in |42j . 



3.1 A family of gauge matrix models 



We will focus on the 2 dimensional version of (2.11). We set 

_ Ai + iA 2 .+ _ Ax- iA 2 , . 

A- -j= , A - -j=—. (3.1) 

Then, one obtains 

S n [A] = j d 2 x({l + Sl 2 )A*A ] *A*A ] + {m 2 -l)A*A*A ] *A ] + 2kA*A ] ). (3.2) 

This action shares some similarities with the 6- vertex modej^] although, as we shall see in 
the following, the entire analysis relies on the choice of a vacuum around which we shall 
perform fluctuations. We will come back to this issue in a while. 

It will be convenient to use the matrix base whose properties relevant for the ensuing 
analysis are recalled in appendix [B| For more details, see |40j . 

The strategy used in this section is standard, akin to the machinery of background 
field method used e.g in [32]: we choose a particular vacuum (the background), expand 
the action around it, fix the background symmetry of the expanded action. The equation 
of motion stemming from (3.2) is 

(3ft 2 - 1)(^U A * A + A * A * A ] ) + 2(1 + Q 2 )A *A ] *A + 2kA = 0. (3.3) 



5 We thank H. Steinacker for this remark. 
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From now on, any solution of (3.3) will be denoted by Z{x). In the following, we will 
consider solutions which respect the symmetries of the classical action (symmetric solutions 
from now on), which were derived in |33j. They have the generic form 



Z„(x) = $i(x z )x^ + $ 2 (x 2 )x M (3.4) 



where <3?i and $2 are suitable functions that have to satisfy Eq. (3.3). Passing to the 
matrix base, one can write equivalently 

Z[x) = -j=- — = Z mn f mn (x) (3-5) 

* m,n£N 



where Z mn can be expressed in terms of $1 and 3>2- This yields 

Z mn = -ia m 5 m+ i jn , Vm, n G N (3.6) 
where the sequence of complex numbers {a m ,m G N} satisfies 

a m [(3fi 2 -l)(|a m+ i| 2 + |o m _i| 2 ) + 2(l + 2 )|a m | 2 + 2 K ] =0, (3.7) 



in view of (3.3). 



This implies, Vm G N, a_i = 0, 

• i) a m = 

• ii) (3Q 2 - l)(|a m+ i| 2 + |a m _i| 2 ) + 2(1 + tt 2 )\a m \ 2 + 2k = 

The first solution corresponds to a trivial vacuum, whereas the second one engenders a 
whole family of symmetric vacua which depend on the range of values of the parameters 
Q and k. 



Setting formally A = Z + (ft, A' = Z* + (jr in (3.2), where (ft can be interpreted as a 
fluctuation around Z, we obtain from the expansion of the Lagrangian, up to an unessential 
constant (we drop from now on the Moyal product symbol *) 

S[<j>,^] = S 2 + S 3 + S 4 (3.8) 

with 

S 2 = y'd 2 x[(/> t (2(l + f] 2 )ZZ t + {m 2 -l)Z ] Z)(ft + 0(2(1 + tf)Z ] Z + {m 2 -l)ZZ ] )(ft ] 

+ (£t((3n 2 - \)zz)$ + <p((3n 2 - \)z^z^)<ft 

+ (m 2 -i)(z^z^4> + z^^Z(/)) + (i + n 2 ){z^z^ + z^z^(j)) + 2K^(j)], (3.9) 



s 3 = J d 2 x [(1 + n 2 ) [z$<ft$ + z^<ft$<ft + z$(ft$ + z^4>U) 

+ (30 2 - \)(Z$$<ft + Z^$<ft$ + Z ] (jxft(ft ] + Z<fxfi<fi] (3.10) 



Si 



y d 2 x [(1 + n 2 )^^ + (3n 2 - i)<ft(ft ] <ft ] 4>] . (3.11) 



Notice that, when expanding around the trivial vacuum, the action (3.8) simplifies con- 
siderably 

S (z=0) [(ft, 0t] = J d 2 x [2as0V + (1 + Q 2 )^^ + (3^ 2 - ( 3 - 12 ) 
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This is the action of a local matrix model, the so-called 6-vertex model[H] (also see 
where the model has been solved). In the following we shall concentrate on the non-trivial 
symmetric vacua which give rise to non-local models. 

The action S[(f>, in Eq. (3.8) is invariant under a background transformation, that 



can be expressed through a nilpotent BRST-like operation, 5 Z , with structure equations 
given by (star Moyal product is understood): 

5 z <p = -i[Z + (f>,C], Sztf = -i[Z ] + 4>\ C], 5 Z Z = 0, 5 Z C = iCC. (3.13) 



The background gauge symmetry can be fixed by adding to (3.8 ) the following gauge-fixing 
action 

Sgf = S z [ d 2 x CQ{4>) = f d 2 x (bG (</>)) - C6 z g{<f>)), (3.14) 



where the structure equations ( |3.13| ) have been supplemented with 5 Z C = b, S z b = 
and G((f>) is the gauge function. The grading and ghost number assignments are as in the 
section 2. 

We find convenient to choose G{4>) = {4> ~ ) so that 

Sqf = J d 2 x {b{4> - <fl) + iC[Z -Z* + <f>-<ft, C}). (3.15) 

In order to simplify slightly the situation, we now integrate over the b field. In the BRST 
language, this amounts to consider a "on-shell" situation (for which now the nilpotency 
of the 5 Z operation is fulfilled modulo the ghost equation of motion together with the 
5^-invariance of the gauge- fixed actiorl 6 |. Doing this generates the constraint <j> = <\y- into 



(3.15) and (3.8). Then, (3.15) becomes 

S GF = j ' d 2 x iC[Z - Z\C] (3.16) 



so that the ghosts decouple as expected. The action (3.8) simplifies to a functional of (f> 
only 

S[</>] = S 2 + S 3 + S 4 (3.17) 

with 

S 2 = J d 2 x [4>{{l + 5ft 2 )(ZZ t + Z^Z) + (30 2 - 1)(ZZ + Z*tf)<t>) 

+2(3ft 2 - l)Z4>Z^<j> + (1 + ^ 2 ){Z<t)Z4> + Z^ZV) + 2K(jxj)] , (3.18) 



S 3 = 80 2 J d 2 x (Z + Z f )0#; S 4 = 4fl 2 J d 2 x <M><$>4>. (3.19) 
Plugging the expansion of <p in the matrix base, <f>(x) = Yl m n ^mnfmni^), into (3.17) gives 



rise to a matrix model of the form 

S[4>] = <t>mn<i>klGmn;kl + Si nt (3.20) 

m,n,fc,ZeN 

where the kinetic operator is given by 



G m n;u = (1 + 5ft 2 ) ^ 5 ml {Z nt z\ k + z\ t Z tk ) + (3ft 2 S m i(Z nt Z tk + Z\ lt z{ 



teN ten 



2(3ft 2 - l)Z ] nk Z lm + (1 + n 2 ){Z nk Z lm + Z\ k z\ m ) + 2nS ml S nk , (3.21) 



6 For a general discussion, see [43] . 
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with Z mn given by (pT6|), Z, 



4 



ia m 5 n +i,m and the cubic and quartic interaction terms 



arc 



Si 



int 



m,p,q,r£N 



Q-r&r+l,m) "I - ^ ^ 4 > mn4 > nk4 > kr4 > rm' 

(3.22) 

Notice that the interaction terms both vanish whenever Q 2 = 0, as it should be with our 
gauge choice, akin to the commutative situation for 2-d QCD in the axial gauge. Note 
also the occurrence of a cubic interaction term. 



The kinetic operator (3.21) can be rewriten as 

Gmn-ki = (1 + bVL 2 )5 m i5 n k{a n a* n+l + a*a ra _i) 

— (3f2 2 — l)(5 m i5 n+ i^-ia n a n+ i + 5mi5 n -i,k+ia* n a n-i — 25 m ^ + i5k+i, n a* n a/) 

— (1 + i} 2 )(5k,n+iSm,i+ia n ai + S n ^k+iSi )m +iana^) + 2K5 m i5 n k (3.23) 



where the constrained by rt3.7h. It can be observed that it involves 2 types of 

terms. The terms proportional to (3f2 2 — 1) are non vanishing when m+n=l+k±2 
while the remaining terms do not vanish whenever m + n = I + k which is similar to 
the so called conservation law for indices that occurs in the matrix formulation of the 
Grosse-Wulkenhaar model |13|. 



3.2 The propagator at the special value Vt 2 = | 

The propagator denoted by P m n-,ki is defined by 

^ G mn;klPlk;sr — &mr$nsi ^ Pnm\lkG kl;r s — $mr$ns' 



(3.24) 



k,l 



k.l 



It depends on the background which corresponds in each case to some solution of (3.7). 

To simplify the ensuing discussion, we now assume that the a m 's are reaQ We also 
assume k ^ 0. The computation can be easily done when the terms proportional to 
(3f2 2 — 1) are absent, i.e when O 2 = ^. The corresponding solution of (3.7), determined 
in [33] is given by: 



1 



-, k < 0, a Tl 



1 



\/-3k, Vm G N. 



Then, (3.23) becomes 

,(1/3) 



G 



mn\kl 



-K)(25 m [5 n k — 5k >n +l5' 



m,l+l 



Jn,k+1 



and satisfies 



G mnll ^ 



m + n = k + I. 



(3.25) 



(3.26) 



(3.27) 



This conservation law acting on the indices implies that Eq. (3.26) depends only on 2 



indices. Indeed, setting n = a — m, k = a — I, with a = m + n = k + l into (3.26) yields 

,(1/3) 



G 



rn,a—m:a—l,l ' 



G 



ml 



\i (25 m i - <W+i - Si m+i), Vm, I G N 



(3.28) 

Gml tO 



where fi 2 := —k. Notice that it does not depend on a. Therefore, we set G 
simplify the notations. 

One observes that G m i is an infinite real symmetric tridiagonal matrix which can be 



related to a Jacobi operator. Therefore, the diagonalization of (3.28) can be achieved by 



therefore setting the arbitrary phases £ m to zero, in the notation of [33] 
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using a suitable family of Jacobi orthogonal polynomials. Note that a similar situation 
arises within the scalar Grosse-Wulkenhaar model [13] as well as in noncommutative scalar 
field theory on 1R? constructed in [21] . It is useful to recall here some technical points that 
will clarify the computation. For more mathematical details, see e.g [46]. 



3.2.1 Jacobi operators 

A Jacobi operator J acting on the Hilbert space £ 2 (N) with canonical orthonormal basis 
{cfcjfceN can be defined as 

(Je) fc = a k e k+1 + b k e k + o fc _ 1 e fc _i, k > 1; 

(Je) = a ei + b e , (3.29) 

where {a k } ke fq&nd {b k } ke fq are infinite sequences of real numbers, with a k > 0, \/k G N. It 
therefore can be represented as an infinite real symmetric tridiagonal matrix. Denoting by 
T>(£ 2 (N)) the dense subset of £ 2 (N) involving all finite linear combinations of the e^'s, one 
can verify that (JX,Y) e 2 = (X,JY) £2 , VX,Y € V(f(N)) where is the usual scalar 
product on ^ 2 (N). Then, J extends to a densely defined symmetric operator on £ 2 (N). 

If in addition J is a bounded operator on V(£ 2 (N)), it extends by continuity to a 
self-adjoint (bounded) operator on 2?(£ 2 (N)). This occurs whenever 

supflafcl) + supfl&fcl) < oo (3.30) 

k k 



which is clearly verified by the operator defined by J m \ = —G m i from (3.28). Then, 
as a corollary of the spectral theorem, the so called Favard theorem [47] (see also [46]) 
guaranties the existence of a set of (real) polynomials {p k (x)} k ^N orthogonal with respect 
to a unique compactly supported measure dfx(x), namely 



(Pm,Pn) ■= / d/j,(x)p m (x)p n (x) = 5 mn , (3.31) 
Jr 

and verifying the following 3-term recurrence relation: 

xp k (x) = a k p k+ i(x) + b k p k (x) + a k ^xp k ^x{x), k>l 

xpo(x) = a pi(x) + b po(x). (3.32) 

Finally, provided some assumptions are verified [47], which will be the case below, the 



above scalar product (3.31) can be expressed as the limit of a particular discrete inner 
product, namely 



N 

(p m ,Pn) = Jim (S^w 2 Ap m {tj)p n (tj)) (3.33) 

where the tj's and wj are respectively called the nodes and the weights. For more details 
see e.g [47]. 

In the following, we will apply this formalism to the diagonalization of the kinetic 



operator in Eq. (3.28) through the determination of the relevant polynomials and related 



measure of integration. 
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3.2.2 Diagonalization of the kinetic operator 



Denoting generically by Afc, k G N the eigenvalues of G mn (3.28), we write it as 



with 



(3.34) 



(3.35) 



where vXnn = T^nm- Then, the combination of Eqs. (|3.26|), (|3.34|) and (|3.35) gives rise to 

(3.36) 



the following 3-term recurrence relation 

Km+l,q + T^m-X,q ~ (2 - Xq^mq = 0, Vm, g G N, 

completed with 



= 0, ^ , g = 1, Vg G N. 

For further convenience, we set 



Pq 



-X q , TZ m (p q ) := TZ mq =, Vg G N. 



Then (3.36) translates into 

U m +l{Pq) + Tl m -l(p q ) = (2 + Pq)Km(Pq), Vm, (/ G N 



(3.37) 
(3.38) 

(3.39) 



with TZ-i(pg) 
equation |^] 



and IZo(pg) = 1 which stems from the following 3-term recurrence 



TZ m+ i(x) + TZ m -i(x) 
ftl(x) 



(2 + ar)7^(x), Vm > 1 
(2 + x)72<,(x) 



(3.40) 



with 72-0 (^) = 1 (JZ-i(x) = 0), evaluated at x = p g , q G N. 

Eq. ( |3.40 1 is of the type given by (3.32). Now, we restrict ourselves to N x N 
submatrices, therefore we impose a cut-off on the indices, namely < m, I, .. . < N — 1 
and define := (—GiL,) in order to make contact with the notation in 3.2.1 Then Eq. 



(3.40) can be cast into the matrix form (matrix product understood) 



J 



N 



( Ko(x) \ ( 

nux) 



+ 



\ 



V n N (x) J 



( K Q {x) \ 



\ n N ^{x) j 



(3.41) 



Eq. ( |3.41 ) readily implies that the eigenvalues of J N are exactly given by the roots of 
TZn(x). Finally, setting 2t = 2 + x in (3.40) yields 



U m+1 (t) + C/ m _i(t) = 2tU m (t), Vm G N*, C7_i = 0, U (x) = 1 



(3.42) 



8 One more intuitive way to obtain the recurrence equation is to compute the secular equation for the 
matrix G m i ( |3.28[ ) when m — 0, m < 1, m < 2,... Calling Dn(x), N = m + 1, the secular equation, 
one observes by induction that Djv+i(a;) + Dn-i(x) = 2xDn{x) and in particular Di = 2 — x, -D2(:r) = 
(2 - i)Di(i) - 1, D 3 (a;) = (2 - a;) Da (a?) - £>i with D = 1. 
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This defines the recurrence equation for the Chebyschev polynomials of 2nd kind 
U m (t) := (m + l) 2 Fi(-m, m + 2; -^), Vm G N, 



(3.43) 



where 2F1 denotes the hypergeometric function. Note that the t7 m (x)'s are a particular 

>«>/^ TT U\ - P^IJA 



family of Jacobi polynomials [Mj P% (x), U m (t) - — ^y- 



Putting all together, we can write 



n m {x) = f{x)U m {^^), Vm 6 N, 



(3.44) 



where the overall function f(x) will be determined in a while so that (3.35) holds true. 

The eigenvalues of and therefore are now entirely determined by the roots of 
U^it). These 9 are given by = cos(^q^), k = 0, 2, N — 1. Then, the eigenvalues 
for the kinetic operator are 



M 2 A^ = 2^ 2 (1 - cos( ^±^)), fc e {0, 2, iV - 1}, 



and satisfy for finite N 
Summarizing our results, we have obtained: 



< fi 2 X k N < 4/x 2 . 



n 



N 
mq 



f(N,q)U m (t") = f(N,q) 



■ r 7r(m+l)( g +l) l 
bm| AT+1 \ 



sin 



77(9+1) -1 

7V+1 J 



0<m,q<N - 1 



(3.45) 
(3.46) 

(3.47) 



where we used C/ m (cosf 



sin((m+l)0) 



To determine the normalization function f(N,q), we first obtain from Eq. (3.35) 

(3.48) 



N-l N-l 

£ KiPS = f(m,N)f(l,N) E U p {tl)U p {tf). 

p=0 p=0 



The sum in the RHS of (3.48) can be computed by adapting the Christoffel-Darboux 
formula [47 1 to the present situation. Indeed, multiplying the recurrence equation (3.32) 
for p m (x) by p m (y) as well as the recurence equation (3.32) for p m (y) by p m (x) and 
summing both, we obtain for any x,y G [—1, 1]. 

N-l 



( x - y) ^2 Pk{x)pk(y) = PN(x)pN-i(y) -PN~i{x)pN{y), x ^y, (3.49) 

k=0 

N-l 

^2iPk(x)) 2 = p' N (x)p N -i(x) -p N _ l {x)p N (x), (3.50) 



fc=0 



where f'(x) denotes the derivative wrt x. Applying these relations to the t7 m 's, (3.49) 
implies that ( |3.48 ) automatically vanishes whenever m / I in view of U]\f(t^) = 0. When 
m = /, we compute 



N-l 



N\ 



N + 1 



(T n+1 (C)Un-i(C)), (3-51) 



((^) 2 " 1) 



p=0 

9 Recall first that, for any ]V G N, a Chebyshev polynomial of 2nd kind E/jv(i) has TV different simple 
roots in [— 1, 11. 
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where the second equality in ( 3.51| ) stems from the relation 

(N + l)T N+1 (x) - xU N {x) 



U' N (x) 



x 2 -l 



and T/v(x) denotes the iV-th order Chebyshev polynomial of first kind 
using T/v(cos#) = cos(N8) in (3.51), we finally arrive at 

N-l 



(3.52) 
By further 



p=0 



sin 



sin 



N+l J 
3r (m+l)7r -i 
N+l 1 



(3.53) 



One easily verifies that the RHS of (3.53) is actually positive for any value of m £ N. 
We finally obtain 



< m = f(N,m)U p (C), 



f(N,m) 



(-l) m (iV + l)- 



sin[^J^] 



sm 



3r (m+l)7T -i 
7V+1 J 



< p,m < N - 1. 



(3.54) 
(3.55) 



Once we have the polynomials which diagonalize the kinetic term we can invert for the 
propagator. Keeping in mind Eqs. (3.24) and (3.28), we set P mn := P m ,a-n;a-i,i where 
a = m + n = k + l. It follows from the above that for fixed N the inverse of G^ n denoted 
by P^ln can be written as 



1 N ~ l 1 



(3.56) 



Taking the limit N — > oo, the comparison of the relation 5 m i = ^2 p T^mpT^i p where the 
7*^ n 's are given by (3.54), (3.55) to the orthogonality relation among the Chebyshev 
polynomials U n 



7T / 

d/j,(x) U m (x)U n (x) = -5 mn , dfi(x) = dxy 1 - x 2 



(3.57) 



permits one to trade the factor f(N,p) 2 (see (3.55)) in P^ n (3.56) for the compactly 
supported integration measure dn{x) ( |3.57 ) in (3.56). 

To conclude this paragraph, we obtain the following rather simple expression for the 
inverse of the kinetic operator (3.28) 

Pmn;kl 
Pml 



1 



J-l 



dx 



l + x 
1 — X 



U m (x)Ui(x). 



(3.58) 



It can be easily checked that (3.24) is verified by combining (|3.58 ) with (3.26) and using 
the orthogonality relation (3.57). Notice that the integral in (3.58) is well-defined leading 
to finite P m i when m and I are finite. 



3.2.3 Computation of the 1-point function 

Finally, let us consider briefly the 1-point function generated by the cubic vertex. The 
computation is standard. Supplementing S[4>] ( |3.20 ) by a source term J^mn^mnJnm, the 
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perturbative expansion can be obtained from the generating functional of the connected 
correlation functions W[J] where 



Z[J] 



n 



S[4']-T. m ,„ 4>mnJ n r, 



,W[J] 



(3.59) 



W[J] = lnZ{0) + W o [J]+ln(l + e- Wo W(e- Sint (sj) _ i) e Wb[J]) ( 3 _ 60) 

Wq[J] = - ^2 JmnPmn;klJkl (3.61) 



m,n,k,l 



and Si n t[(j)] collects the interaction terms of the action, which can be read off from Eq. 
(3.22), evaluated at ft 2 = §. The ] propagator P m n-kl is given by Eq. (13.581) and the Jnm's 



are sources. Expanding the logarithm as a formal series yields all the connected diagrams 
while the effective action T[(f>] is obtained from W[J] by Legendre transform 



'mn' J nm 



W[J], 



SW[J] , 

5J n 



J=0- 



(3.62) 



From the formal expansion of (13.60) we obtain 



W^I^I — ^ v mn((Plm;kl + Pkl;lm)(Pnk;cdJcd + JgbPgb; 
+ {Plm;nk + Pnk;lm)(Pkl;cdJcd + JabPab;kl) 
~\~ (Pkl;nk ~i~ Pnk;kl){.Plm;cdJcd J abPab;lm)) i 



nk J 



with a 



3fi 2 and 



Jm+l,n 



3m,n- 



1- 



(3.63) 



(3.64) 



Combining (3.63) with the relevant term in the expansion solving the 2nd equation in 

Jmn = G nrn -kl4>kl + •••• (3.65) 



(3.62 ) given by 



k,l 



where G nm -ki is given by (|3.26|), we obtain 



a 
2 



ab:nk J 



^ ^ ^mn y^{Plm;kl Pkl;lm) {Pnk;cdG dc;rs ~i~ Gl} a - rs P a , 
~i~ (Plm;nk Pnk;lm)i.Pkl;cdGdc;rs Cftajrs-fafcjfci) 

+ (Pkl;nk + Pnk;kl)(Plm;cdGdc;rs + Gb a -rsPab;lm) 4>rs- (3.66) 



Now combining (3.66) with (3.58) and 



nl)9ml 



{Plm;cdGdc;rs ~\~ Gl )a - rs P a f,-lrn)4 > rs — 20 m ; 

we can write T 1 [4>] into the form 

rl [0] = a V mn(5mk(Pll + Pkm)<Pkn + S n l(Pkk + Pn 

m,n,k,l 

,n+kPlk + dn+k,rn+lPnin)4>lk) ■ 

By using the Kronecker delta symbols, r 1 [</>] can be cast into the form 

r 1 ^] = <y /X^Pu - Pi,i+i + Pkk + Pk+i,k+i - Pk,k+i)(4>k,k+i - 4>k+i,k 



(3.67) 



(3.68) 
(3.69) 
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It is divergent, with typical divergence 



rjU+i~l>i^-A!+i) 



(3.70) 



This signals likely that the vacuum (3.25) is not stable against quantum fluctuations. 



Indeed, from the structure of r 1 [^>] as in Eq. (3.69), it seems difficult (if possible at all) 
to absorb the divergences by the natural set of counterterms depending on two arbitrary 
parameters that are generated from the action (3.2) (when Q 2 = |) and further expanded 



around the vacuum, akin to the linear-sigma model. 



4 Discussion and conclusion 

A lot of information can be extracted from the propagator whose fate, and consequently the 
fate of the corresponding matrix model, is completely determined by the chosen vacuum. 
The UV and IR region can be identified from the spectrum of the kinetic operator. First, 



by taking the limit N — > oo of the spectrum given in (3.45), one observes that 



lim n 2 \to = 0, hm f^X^N-l = V- (4.1) 

iV— >oo N-^oo 

Then, from an overall rescaling of the initial action by a factor g~ 2 with mass dimension 
[g] = 1 so that [/i 2 ] = 2 (and [A] = 1), it is natural to define the UV region as the one 
corresponding to large indices while the IR domain corresponds naturally to low indices, 
say, m = 0, 1. 

Next, we observe that the operator on £ 2 (N) defined by the matrix elements P m i leads 
to a propagator which does not decay at large separation \m — l\. This can be already 
realized by expressing the matrix model in the field variables diagonalizing the kinetic 



operator, i.e the propagation base in the physics language. From (3.45) and taking the 
large N limit, one readily finds that the eigenvalues of the propagator // _2 A£° increase 
as the index k increases which can be interpreted as signaling a lack of UV decay for 
the propagator. This can be verified directly on (3.58) by computing numerically P m i. 



Unsuppressed correlations P m i may exist at large separation \m— 1\. This is to be compared 
with the situation for the Grosse-Wulkenhaar model for which there is a suppression of 
the correlations at large separation which then mimics a kind a quasi-locality, and was one 
of the main ingredients leading to the renormalizability of the model. Such a behaviour 
indicates that the matrix model under consideration, at the value f2 = | is highly non 
local. 

We note that the interpretation of the action for the matrix model as the spectral 
action related to a finite volume spectral triple is not obvious, if possible at all. Such 
a triple, as introduced in [29], is characterized by a self-adjoint Dirac operator D with 
compact resolvent operator whereas the defining algebra of the triple is not unital. In the 
present case, one should have D 2 = G (G being the kinetic operator) from the spectral 
action computation. But G is bounded, so that G + I is also bounded. But then, (G + I) -1 
is not compact. Hence, the resolvent operator is not compact. 

The lesson to be drawn is that requiring the propagator to decay at large indices, i.e 
in the UV, singles out one specific family of vacua among those symmetric vacua classified 
in [33]. Let us analyze this point in some detail. In |33j the symmetric vacua have been 
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classified accordinj 


I to 








o < n 2 


1 

<3 ' 










<1 , 


u 2 - 




m ), k < 0, r < -1 


SI 2 


= 1 , 


U 3 m = 


-fd-(-l) 


~ m ), K < 



with 



r 



i + n 2 + V80 2 (i -^ 2 ) 



(4.2) 

(4.3) 
(4.4) 

(4.5) 



1-3^2 

/u^ and t the arbitrary phase set to zero as in the previous section), besides the 



solution (3.25) for Q 



i 



that we have analyzed in detail, which leads to a propagator 



with an unsuitable decay. A similar conclusion applies to the solutions (4.3) and (4.4) 
Indeed, one observes that lim m _ s . 00 |ii^J = and lim-m^oo \u^ n \ < In view of (3.23) 
this corresponds to a bounded self-adjoint operator. Thus, the real spectrum, Spec(G) 
is a compact subset of [—UGH, ||G||]. This should not give rise to a decaying propagator 



This observation does not apply to (4.2) for which the kinetic operator is unbounded, in 
view of 

I — i — i m 

V n m = a m ~ r 2 , m — > 00. (4-6) 
Then, the corresponding propagator is likely to have a suitable decay behaviour. This 



altogether singles out the family of symmetric vacua defined by (4.2). However, Eq. (4.6) 



indicates that the vacuum defined by Eq. (4.2) does not belong to the Moyal algebra, i.e 



the multiplier algebra recalled in subsection 2.1 This would exclude the vacuum solution 
Eq. (4.2) if one insists to preserve the present definition ofR 2 ,. Thus, one should conclude 
that the interpretation of Moyal gauge theory as a matrix model is problematic, at least 
when using symmetric vacua. However, enlarging to a new set including the solution 



(4.2) provides a way to overcome this obstruction. We think that this is an interesting 



possibility that deserves further investigations. 



We notice that there is a formal similarity between the classical action (3.2) and the 



action describing the 6- vertex model [2]. However, expanding (3.2) around a non-zero 
vacuum gives rise to a different model (the 6- vertex model corresponds to a zero vacuum). 
Once the vacuum (background) is chosen, here among the symmetric non-trivial vacua of 
the classical action, the kinetic part is entirely determined, up to gauge fixing, by the part 
of the expanded action which is quadratic in the fluctuations. The kinetic part therefore 
depends on the vacuum solution chosen, e.g. 



(3.6). 



We have studied a class of noncommutative gauge theories elaborated in [26] on 2-d 
Moyal space from the viewpoint of matrix models. We have explored some related prop- 
erties beyond the classical order. Expanding the action around symmetric vacua classified 
in [33] generates non local matrix models with polynomial interaction terms. For a par- 
ticular symmetric vacuum, we have shown that the kinetic operator is a Jacobi operator. 
The computation of the propagator has been carried out. The resulting propagator can be 
expressed in terms of Chebyschev polynomials of second kind. We have shown that non 
vanishing correlations exist at large separations. For such particular vacuum a divergent 
1-point function appears that seems difficult to absorb by the set of symmetric countert- 
erms, signaling possibly a loss of symmetry at the quantum level. The quantum stability 
of the vacuum is briefly discussed. 

From general spectral considerations on the kinetic operators stemming from the other 
classes of vacua determined in [33], we have singled out a particular class of symmetric 
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vacua that should lead to fast decaying propagators at large separation corresponding 
to quasi-local matrix models. The latter class of vacua is an interesting possibility to 
obtain consistent (quantum) matrix models from the present scheme and deserves further 
investigations. In particular, one should determine whether or not these vacuum solutions 
can be reconciled reasonably with algebraic structures related to Moyal noncommutative 
geometry and compute the propagator from a non Jacobi kinetic operator. 
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on recent developments in matrix models. Discussions with D. Blaschke and H. Grosse 
at various stages of this work are gratefully acknowledged. One of us (JCW) thanks H. 
Grosse for having pointed out past works related to spectral properties of operators. 

A Vertex functions at Q = 

It is convenient to define p A k = p^Q^uku and p^ = Q^py. We list below the vertex 
functions p*| for the Yang-Mills theory M 2 , (Q = 0) in the Landau gauge. The propagator 
for the is G^ v {p) = 6^/p 2 . The ghost propagator is G cc (p) = 1/p 2 . 

k\ A k 2 

Vapjifa, k 2 , k s ) = -z2sin( — - — ) [(k 2 - fci) 7 <5 a /3 + (h - h)pd ai + (k 3 - k 2 ) a Sp 1 ] , (A.l) 
Vap 1 6( k i,k 2 ,k 3 ,k 4 ;) = -4 [(SarySpg - 5 aS 5p y ) sin( fcl ^ ^ 2 ) s i n ( k?J ^ ki ) 



, rx x x x \ ■ ,h /\k4 x . ,k 2 /\k 3 . k 3 A ki , . M A fc 4 -, 
+ [p a p5 lS ~ OayOps) sm( — ) sm( — ) + {d aS 6p 1 - d a pd lS ) sm( — ) sm( )J , 

(A.3) 

V cc „(k l ,k 2 ,k 3 ) =i2A; 1/1 sin(^^). (A.4) 

The IR behaviour of the correlation functions can be conveniently controled by using the 
integrals given in |27j : 

/d D k e ik P 
(2ir) D (fc 2 + m 2 ) N = a N,DM N _n(m\p\), (A.5) 



/dPk k k e ikp 
(2^p (k 2 +m 2 ) N = aN M S ^ M N-l-fi m \p\) ~ P»PvM n _ 2 _d (m\p\)), 

(A.6) 

where 

r(iV)vrf ; (m 2 f 



a N,D = n ; M Q (m\p\) = - r -^{m\p\) Q K Q {m\p\) (A.7) 



10 Momentum conservation is understood. All the momenta are incoming 
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in which Kq(z) is the modified second kind Bessel function of order Q £ Z. Recall that 
one has 

K_ Q (z) = K Q {z); lim z v K v (z) = 2 v ~ 1 Y{v) i v > (A.8) 
so that the following asymptotic expansion holds true: 

M- Q (m\p\)^2Q- lT ^, Q>0. (A.9) 

B The matrix base 

The matrix base used in the section [3] is the family of Wigner transition eigenfunctions of 
the harmonic oscillator [50] {/ m n(^)}m,neN C 5 which can be expressed as 

fmn = 1 * f00 * (B.l) 

(mln\6 m+n )2 

with /oo = 2e~~s zz and a* n := a-ka-ka* ...*a (n times). The following relations hold true 
{fmn * f Pq ){x) = S np f mq (x), / d 2 xf mn (x) = 2ir95 mn , f mn (x) = f nm {x), Vm, n,p,g£N 



(B.2) 

Any element of the Mg can be writen as a = Ylm n a mnfmn so that the ^-product between 
2 elements of Mg is mapped into a product of 2 (infinite) matrices, namely (a-kb)(x) = 
J2 m ,n c mnfmn(x), c mn = J2 p a mpbpn- Further usefull formulas are 

/oo = n^^- 1 ) * /oo, /oo * ^ n * z = n9f oo * z< n ~^ ,n>l (B.3) 
while both RHS are zero for n = 0. 
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